Mimicking interacting relativistic theories with stationary pulses of light 
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One of the most well known relativistic field theory models is the Thirring model (TM). Its real- 
ization can demonstrate the famous prediction for the renormalization of mass due to interactions. 
However, experimental verification of the latter requires complex accelerator experiments whereas 
analytical solutions of the model can be extremely cumbersome to obtain. In this work, following 
Feynman's original proposal, we propose a alternative quantum system as a simulator of the TM 
dynamics. Here the relativistic particles are mimicked, counter-intuitively, by polarized photons 
in a quantum nonlinear medium. We show that the entire set of regimes of the Thirring model - 
bosonic or fermionic, and massless or massive - can be faithfully reproduced using coherent light 
trapping techniques. The sought after correlations' scalings can be extracted by simple probing of 
the coherence functions of the light using standard optical techniques. 



Introduction 

The Thirring model is the simplest relativistic quan- 
tum field theory that can describe the self interaction of 
a Dirac field in 1+1 dimensions pQ. Various techniques 
have been developed to predict its correlations function 
behaviour with conclusive results only for the massless 
fermonic case [2j [3] . For the opposite case of fermions 
with mass, approximations involving the form factors ap- 
proach have been developed for the n-point correlations 
whereas a full analytical solution still evades the efforts 
of seminal theoretical physicists [3H8]. The theory of the 
bosonic regime and the correlation function behavior is 
less developed with efforts mainly concentrated in the 
solving the classical case [9| [TO]. 

The possibility to physically simulate the Thirring 
model dynamics in a different controllable quantum sys- 
tem is thus extremely interesting. An ideal quantum sim- 
ulation should be capable of reproducing on demand the 
bosonic or the fermionic cases with tunable interactions 
and mass terms. In addition direct access to the corre- 
lation function behaviour will be an advantage. Such a 
system would not only shed light in the unknown parts of 
the Thirring model itself, but motivate works for prob- 
ing of the behaviour of other relativistic quantum field 
theories which are beyond the realm of numerical and 
experimental techniques. In this work we describe ex- 
actly such a setup where the TM is generated counter- 
intuitively in nonlinear optical system. Here polarized 
photons are mimicking the behaviour of the relativis- 
tic fermions where interactions and the mass terms are 
generated and controlled using quantum slow light tech- 
niques. 

Recent works in quantum simulations of relativistic 
theories focused on free fermion simulations using cold 



atoms and ion systems [TTHT3] . An interacting fields 
proposal is described in a recent seminal work where 
fermions on an optical lattice are used [14] and more 
recently a connection between a variational simulation of 
fields to a generic cavity QED setup was made [15 j. The 
present work differs in many aspects to all previous quan- 
tum simulations. First and foremost it is not based on 
cold atoms or particles with mass in general, but on a 
fundamentally different and somewhat counter intuitive 
approach using photons in a nonlinear optical medium. 
In addition, our optical setup being a natural contin- 
uum system, has the ability to reproduce the continuous 
Thirring model directly without resorting to a lattice ap- 
proach. Finally, both the bosonic and the fermionic cases 
can be generated by controlling the relevant optical inter- 
actions, and in both massless and massive regimes. In our 
proposal the relativistic dynamics (linear dispersion) and 
the necessary strong interactions are generated and con- 
trolled by employing electromagnetically induced trans- 
parency (EIT) techniques [16HT8] .The scalings of the cor- 
relation functions of the TM, for any regime of interac- 
tions, could be simply probed by analyzing the quantum 
optical coherence functions of the outgoing photons as 
they exit the medium. 

We start our analysis by showing the possibility to 
generate a nonlinear Dirac type of Hamiltonian using 
stationary polarized pulses of light. We then analyze 
the parameter regime to tune the system to a hardcore 
regime reproducing the FTM dynamics and continue by 
showing that both the massless and massive cases are 
reproducible. We conclude by discussing how one could 
probe the scaling of correlation functions for any regime 
of interactions and mass values using standard quantum 
optical measurements. In this last part, we calculate ana- 
lytically for demonstration purposes the correlation func- 
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tion behaviour for the case of the massless fermionic TM 
where an exact solution is known. This could be used to 
gauge or first test our simulator before it is used in the 
territory of theories without known solutions. 



Photons as interacting Dirac fermions: The system 

Our proposal is based on exploiting the available huge 
photonic nonlinearities that can be generated in specific 
quantum optical setups. More specifically, we envisage 
the use of a highly nonlinear waveguide where the nec- 
essary nonlinearity will emerge through the strong inter- 
action of the propagating photons to existing emitters in 
the waveguide. Recent efforts have developed two simi- 
lar setups in this direction, both capable of implementing 
our proposal with either current or near future platforms. 
In these experiments, cold atomic ensembles are brought 
close to the surface of a tapered fiber [19 or are loaded 
inside the core of a hollow-core waveguide {2C[ as shown 
in Fig. [l](a). The available huge photon nonlinearities 
due to the (EIT) effect can be used to create situations 
where the trapped photons obey TM dynamics. 

We start by assuming two probe quantum fields of few 
photons each, labeled as E s , E s >, are tuned to propagate 
in the medium. s, s' =t, I will denote their different po- 
larizations ( later to be representing the spins of the sim- 
ulated relativistic particles). The waveguide can also be 
illuminated by two pairs of counter-propagating control 
lasers with Rabi- frequencies Q s ,± with ± denoting right- 
and left-moving directions. The atoms placed inside or 
close to the surface of the waveguide (depends on the im- 
plementation) have a ground state \a) and excited states 
|6, s) and \d,s,s f ) Each of the probe pulses is paired with 
an appropriate classical field and is tuned to an atomic 
transition such that a typical A (EIT) configuration is 
setup (see Fig. 1(c) and Methods). 

The process to generate the TM dynamics in our pho- 
tonic system is as follows: First, the laser-cooled atoms 
are moved into position so they will interact strongly with 
incident quantum light fields as done in p~9j [20] . Initially, 
resonant to the corresponding transitions, the classical 
optical pulses with opposite polarizations are sent in from 
one direction, say the left side. They are injected into the 
waveguide with the co-propagating classical control fields 
f^ 5+ (£) and Q^+(t) initially turned on. As soon as the 
two quantum pulses completely enter into the waveguide, 
the classical fields are adiabatically turned off, converting 
the quantum probes into coherent atomic excitations as 
done routinely in slow-light experiments. We then adia- 
batically switch on both fl s ,+ and £l s ,_ from two sides. 
The probe pulses become trapped due to the effective 
Bragg scattering from the stationary classical waves as 
analyzed in [18 , 21, 22 . At this stage the pulses are non- 
interacting with the photons expanding freely due to the 
dispersion. By slowly shifting the ^-levels we show that 




FIG. 1: An illustration of the possible realization of our pro- 
posal. The quantum pulses E\,+ and enter the non- 
linear waveguide with different polarizations where they are 
trapped and made to strongly interact. Q^,± and £l±,± are 
two pairs of classical lasers helping to trap and control the 
pulses. In this regime, the polarized pulses can be made 
to mimick interacting fermions of opposite spins following 
the TM dynamics. The scheme is based in engineering the 
photons' interaction with atoms in the waveguide using EIT 
techniques. The appropriate relativistic dispersion relation 
and the relevant TM mass and interaction terms, are fully 
tunable by tuning optical parameters such as atom-photon 
detunings and laser strengths. Once the strongly interacting 
regime is reached, the pulses can be released with all corre- 
lations established during the interaction period coherently 
transfered over to outgoing photon pulses. The latter can be 
measured using standard quantum optics techniques. At the 
insets the atomic level structure and a possible implementa- 
tion using hollow-core fiber or tapered fibers interacting with 
cold atomic gases. 



the effective masses of the excitations can be kept con- 
stant whereas the effective intra- and interspecies repul- 
sions are increased. This drives the system into a strongly 
interacting regime. This dynamic evolution is possible by 
keeping for example the corresponding photon detunings 
A s constant while shifting the <i-level. Once this cor- 
related state is achieved, the counter-propagating con- 
trol fields fl s ,- responsible for trapping the polaritons 
are slowly turned off. This will release the corresponding 
quasi-particles by turning them to propagating photons 
which will then exit the fiber. As all correlations estab- 
lished in the previous step are retained, the propagating 
wave packets will carry the correlation of the TM dynam- 
ics taking place in the medium outside, where they can 
be probed by appropriately placed photon detectors. 

A Thirring model of polarized photons: In the 

trapped regime the dark-state polaritons emerge in the 
form \I/ S ,+ = jy^- E s ^+ and ^ s ,- = Es,-i which are 
also propagating in both directions (see Methods). The 
symmetric combination \I/ S = a s?+ \I/ S5+ + a S5 _\I/ S) _ with 
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a Sj ± = ^2 — — forms a stationary wave obeying the 
following tunable nonlinear dynamics (see Methods) 

ihd t ^ s = -- d 2 z y s + ihr) s d z y s + hn y-s 

+Xss^ 2 s + x i5 *|*i* 8 + noise (1) 

where the optically tunable inter- and intra- species in- 
teractions are given by 

Xss = 4tef s /(A ss n z ) (2) 
Xs§ =2m 2 a [2 + co8(tp i -ip a )]/(A ai n z ) (3) 

The mass of the particles in the relativistic regime is 
mo, s = — M^o/^and r] s — —2v s cos2(f s . We note here 
that the above mass is different from the "dressed" mass 
due to the interactions. 

The quadratic dispersion term corresponding to the 
"classical" kinetic energy can be ignored in the regime of 
interactions we operate (Methods). In this case, if we de- 
fine a spinor field \l/ = (\f^ , ^|) T , the Hamiltonian gen- 
erating the last equation becomes the bosonic Thirring 
model 

H = J dz[-ih\rj\'Y 1 d z + m r? 2 )^ 

+ ^^^^^s + |*7 /x **7/x*] 5 ( 4 ) 

s 

where 4?= ^70 and the spinor fields satisfy the commu- 
tation relations [W a (z) , * a / (z')\ = 0, (z) , \E^, (z f )] = 
5(z — z')5 ss > . The gamma matrices 7^ for \i = 0, 1 are cho- 
sen as 70 = cr Xl 71 = i<jy with covariant gamma matrices 
7 = 70, 7 1 = —71, and an additional matrix 75 = 7071 • 

A fermionic TM with tunable mass: The Hamil- 
tonian in Eq.(4) describes bosonic particles made of sta- 
tionary light-matter excitations. In order to realize the 
FTM we will tune to the hardcore regime where the same 
species interaction is much larger than all other interac- 
tion terms. In Fig. 2a we calculate and plot the ratio 
of their strength as a function of the controllable sin- 
gle photon detunings. We see that by appropriate tun- 
ing from the corresponding upper atomic states ( tun- 
ing the laser towards or away from the corresponding 
states) the system enters a hardcore regime where the 
same species interaction is much larger than any other 
term in the Hamiltonian. The particles will then be- 
have as effective fermions in all aspects-identical density 
density correlations [23, 24 and spectra but different first 
order correlations due to the sign issue. We note here 
that as we will discuss later, our optical simulator will be 
naturally probing intensity-intensity correlations which 
nicely fits with the above point. 

For this regime, we also calculate and plot the individ- 
ual interaction strengths, denoted Pl s (/?*-), as a function 
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FIG. 2: Fig 2a: The ratio of inter- to intra species repulsion 
as a function of the relevant single photon detunings and the 
relevant "bosonic" and "fermionic" regimes; Fig 2b(c): The 
ratio of intra(inter) species interaction to kinetic energy, de- 
noted by /3l s (Pig), as a function of the corresponding single 
photon detuning 

of the kinetic energy of the particles for the same values 
of the tunable optical detunings (Fig. 2b and 2c). We 
see that when the system is optically tuned to the hard- 
core regime, the intra-species repulsion is much larger not 
only the inter-species interaction (Fig2 a) but the kinetic 
energy as well (Fig2b,c). In plotting (Fig2 a) we have 
set mo = with ^0 turned off. We note here the simple 
way of tuning the rest mass of the particles by simply 
controlling the "connecting" laser. 

Probing the TM correlation scalings with photon 
detectors 

In the following we describe how one could probe the 
scaling behaviour of the correlation functions and from 
that infer the model's properties and the fascinating 
renormalization of mass due to the interactions. We use 
the case of the massless fermionic TM as an example to 
illustrate our approach. Of course in this case the spec- 
trum of the model and the scattering matrix was explic- 
itly evaluated by the Bethe Ansatz [4] and is known. Nev- 
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ertheless, this regime serves as an excellent check for our 
quantum simulator. Ultimately, one would like to move 
to the relatively unexplored regimes of massive or bosonic 
TM as well as other interesting relativistic theories [25] 

In the massless FTM the long-scaled n-point correla- 
tion as a function of the distance z — z' and interac- 
tion strength x/ \v\ was calculated by Coleman in [2] and 
reads: 



(oin*i(^)^(^)*i(^)*T^)io) 



l\ 2 "A 2 rU- (z-Zjfiz'-z'/M* 



[i+^r 1 



(5) 



One can extract the two-point correlation with (n — 1): 



<0|*} (z 1 )* l (^)*| (^)*t (*i)|0> 
4 



A 2 (z 1 -z[ Y 



rr\ V \l 



(6) 



A is the momentum cutoff, which we can estimate for our 
finite size waveguide of lentgh L containing N interact- 
ing particles, to be A = 7rn p h sin (xl M) / (x/ \ v\) •> where 
n ph — N/L [4]. In Fig. 3(a) we plot the dependance of 
the cutoff on the ratio of interactions to kinetic energy 
for the regime < xl M < ^ where modes with un- 
physically high energy can be excluded. This regime of 
interactions corresponds to our optical simulator in sin- 
gle photon detuning A ss / ranging up to a few T (which 
is within the required conditions for our slow light setup 
dynamics [El [17]). In Fig. 3(b), we numerically calculate 
and plot the two-point correlation as expressed in Eq. (|6| 
as a function of the distance normalized to the inverse of 
the photonic density for maximum interactions. In these 
plots, as also in the previous one in Fig. 2, we have as- 
sumed the typical values of the optical parameters for 
the slow light setup, ie. |cos2<^ s | = 0.004, n z = 10 7 m _1 , 
= 0.2r, and U s ~ 1.5I\ < xl \v\ < * 
To detect the above correlation in our optical 
proposal, is enough to observe that it is in fact 
an effective transverse spin-spin correlation function 
(S+ fa) S~ &)) with S+Ozi) = tf}^!)*^) and 

$~ ( z i) ~ ^| ( z i) ( z i)- I n our case ? one would re- 
lease and map the polaritons to photons by turning off 
and then using polarization wave plates to convert 
\I/ S into their superpositions ^ x ,± and ^ y ,± as the pulses 
exit the medium which then can be measured by probing 
intensity-intensity correlations (see Methods). The lat- 
ter could be done using standard Handbury-Brown and 
Twiss type of techniques for different times as time cor- 
relations for the propagating photons here map to space 
of the interacting particles. 
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FIG. 3: The momentum cutoff A in unit of 7m p h (a) and log- 
scaled two-point correlations (b) of the massive FTM. As the 
distances z — z in (b) increases, the correlations decrease. We 
simulate the system that each initial quantum pulse contains 
10 photons and spreads over after completely entering into 
lcm-length fiber. 



Conclusions: We have shown how to physically sim- 
ulate the Thirring model dynamics in a photonic setup 
where fermions are mimicked by stationary polarized 
pulses of light in a strongly nonlinear EIT medium. We 
showed in detail how using coherent light trapping tech- 
niques, the photons' dispersion relation can be tuned to 
the relativistic regime and how to create optically tun- 
able interactions and mass terms. The capability to nat- 
urally probe the scaling behaviour of the TM correlation 
functions using standard optical measurements on the 
emitted pulses, makes our approach complementary to 
alternative simulations based in cold atom setups. The 
possibility to first test the device in regimes where an- 
alytical solutions are known (massless regime) and then 
probe regimes with unknown solution (massive FTM and 
bosonic TM) or other interacting relativistic QFTs makes 
this extremely interesting in our opinion. Finally, the 
open nature of our setup and the interaction with the 
environment through photon loss and driving, should al- 
low for the possibility to investigate QFT theories at fi- 
nite temperature, an extremely fascinating and well un- 
explored area of research. We hope that the above char- 
acteristics combined with the simplicity of our approach 
will inspire further experimental and theoretical investi- 
gations in using photonic systems for quantum simula- 
tions of exotic theories. 
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METHODS 



Nonlinear dynamics of relativistic stationary 
polaritons 



The Hamiltonian of the system reads H 
-h J n z dz(hA + + he), where 



?' &d,s,s' :d,s,. 



(7) 



= ^2(g S Vb jS ;a + ^2gss>0-d,s,s>;c,s)Es + h.C, (8) 
s s' 

= ^CT C ,s;6,s^s + ^0^c,t;c,| + h.C.. (9) 
s 

Here, the quantum filed is composed of two counter- 



propagating components, E s 



*Q,sZ _|_ 



E s -(z,t)e lk Q' sZ . Similarly, control fields Q s = 
n s ]+(z,t)e ikc > sZ + ^ s ,_(z,t)e- ifcc ' s2 , where fi Sj ± are 
slowly varying parameters of z and t. kQ jS and /cc, s de- 
note the wave vectors corresponding to central frequen- 
cies ujq jS and co>c,s of E Si ± and fi Sj ±, respectively The 
collective and continuous operators o^ v = a fl . u (z : t) de- 
scribe the averages of the flip operators (v\ over atoms 
in a small region around z. n z is the atomic density 
g s and g ss > are coupling strengths between the quantum 
fields and atoms, while A s and A ss > are one-photon de- 
tunings for corresponding transitions. For simplicity, we 
assume that g s = g ss > = g. 

As done usually in slow light works [T7| H8] . we con- 
sider that atoms are initialized in the ground state \a). 
In the trapped regime the dark-state polaritons emerge 



in the form \I/ S ,+ = 



E s ^+ and \I/ S 



which are also propagating in both directions. Here, 
we have considered the slow-light limit where the con- 
trol fields are tuned to be very weak and the polari- 
tons are mostly spin excitations. The imbalance be- 
tween two counter-propogating components of control 

fields is described by tan 2 ip s = _ / + . Another 

2 

angular parameter is introduced as tan 2 S = g 2 n z /Q s 
with Ti 2 s = (^ 2 + + Vt 2 s _)/2. The reduced light veloc- 
ity in the nonlinear medium is v s = v/(7rtan 2 9 S ), where 
v = co>q 5S //cq 5S is the light speed in the empty medium. 
The dark-state polaritons dynamics is determined by the 
Maxwell-Bloch equations of quantum field 

(d t ± vd z )E s ,± = in z g(a a]bi8] ± + ^ cr c ,s;d, s ,s' ; ±) (10) 



Substituting the dark-state operators \£ s? ± and adi- 
abatically eliminating the fast decaying atomic opera- 
tors and ignoring high-frequency oscillation terms, the 
Maxwell-Bloch equations become 



d t ^s + vd z (a 3i +V a - a s ^ s + 2a 3i +ot a -A a ) 

tan 2 S ^ T .tan 2 6L . 2q 2 T + T 
+ z— ^fio** - z^-***,*, 

a 2 

-z^-[2 + cos(^-^)] 

dtA s + vd*(2tf s + a s -A s - a s ,+A s ) 

A A s A - s 



(ii) 



Here 5,5 =t>i and s ^ s. We keep the symmetric 
combination \I/ S = a S)+ \I/ S5+ + a S) _\I/ S) _ with a S5 ± = 

^2 — ^2 — while the antisymmetric combination A s pro- 
portional to — ^s,- can De adiabatically eliminated 
in the limit of large optical depth since the pulse match- 
ing phenomenon gives — ^ s ,- —> [18]. 

The evolution equations for dark-state polaritons then 
become 



h 2 



2rn r] 



-d 2 y s +ihri s d z y s + hnoy-s 



+X88*t*s + X88*l*8*8 + noise (12) 

with interspecies interaction strength given by 
2 

Xss = ^hQ s /(A ss n z ), while the intraspecies interaction 
strength is x s s = 2h ^s \ 2 + cos (v* ~ Vs)] /{A 8 - 8 n z ). 
m nr;S = — Ml s /[4sin 2 (2(/? s )v 2 A s ] is the polariton mass. 
The mass of the particles in the relativistic regime 
is mo, s = —fc^o/risi which is different from the non- 
relativistic polariton mass m nT ^ Sl and r] s = —2v s cos2(p s . 
Note that the above is different from the "dressed" mass 
due to the interactions. 



The bosonic Thirring model 

We will now show that we can tune our system such 
that the first term on the r.h.s. of the equation, the 
quadratic kinetic energy term, to be much smaller than 
the linear relativistic one (second term). For this we first 
set 77^ = — 77^ = r] to simplify the expressions. The latter 
is possible by tuning 0^ = 0^ through the slow light ve- 
locities for each species = v± and cos 2^ = — cos 2^. 
These correspond to different ratio Rabi frequencies of 
the classical lasers ft 2 + and ^ 2 _ set by tan 2 cp s . The 
non-relativistic kinetic energy reads 

£ nr , s = (Psf /(2m nr , s ) = 2ftsin 2 (2(^K 2 |A S | /[U] (z s f) 

(13) 

and corresponds to the first quadratic dispersion rela- 
tion term above. The relativistic one (second term in 
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r.h.s. of Eq. 6) is, E T)S = \rj s \P s , where P s is the charac- 
teristic momentum of the polaritons. One can estimate 
P s ~ h/z s based on the typical spatial extent of the 
polaritons z 3 , which can range between L for linear or 
weakly interacting photons to L/N p h jS in the Tonks gas 
regime. The latter ranges from L to L/N p ^ s ( weak or 
linear to Tonks gas regime) with 7V p h, s is the s-type pho- 
ton number. Setting their ratio to we can easily see 
that the latter is simple tunable function of one-photon 
detuning: 



ft 



Enr,s _ sm 2 (2(p 8 )v 3 \A S 

^r,s |cQS 2(p 3 \ Z s Q 2 g 



(14) 



Fixing the mixing angle |cos2(/? s | = 0.004, the atomic 
and photonic densities n p h, s = ^ P h = 10 3 m _1 , the coop- 
erativity factor Tid = 0.2r, v s — lOOm/s and Q s ~ 1.5T, 
Pk,s'-n we see that the ratio can be as small as 5% for sim- 
ilar values of photon detuning A s /r from 0.05 to 0.08, 
which allows us to neglect the contribution of E nr ^ s in 



the Hamiltonian Eq. (15). The nonlinear equation then 



becomes the bosonic TM: 



H = [ dz[-ih\ri\jid z + m r/ 2 )^ 



•Xl^*^*»*» + f*^**^*]. ( 15 ) 



where ^=^70 and \l> = (^ t ,^) T with the 
spinor fields satisfying the commutation relations 
[*.(*),*✓(*')] = 0, [* s (*),**,(*/)] = 5{z-z')5 ss ,. 
The gamma matrices 7 M for p = 0, 1 are chosen as 

0": 

7o, 7 

We have also set mo = mo,| = mo,;, which is achiev- 
able by tuning 77^ = —77;. Furthermore we can also have 
m nr , t m nr ,; m nr and 2\n = 2 Xn = X by appropri- 
ately tuning the single and two photon detunings. 



7o 
7° 



x , 71 = i&y with covariant gamma matrices 
1 — —71, and an additional matrix 75 = 7071 • 



Probing of the TM correlations optically 

The two-point correlation function of massless FTM 
is an effective transverse spin-spin correlation function 
(S+ (z ± ) S~ (z[)) with = and 

S~ (z[) = ^{ (z[) ^ (z[). We could probe this by mea- 
suring correlations between the two polaritons species 
densities p x ,± = ^ x ,±^x,± an d p y ,± = ^ y ,±^y,± with 
tfa.i = (*t ± *^)/V2 and * yi± = (* t T i*0/>/2. 
We can express S + = ^j^; = [(Ae,+ ~ Px-) + 

KPv,+ ~ Py,-)}/ 2 and S~ = = - Px,-) ~ 

i(p y ,+ — p y -)}/2, which indicates the correlation func- 
tion (S + (zi) S~ (z[)) of massless FTM as a combination 
of density-density correlations (p a ,± ( z i) Pa',± ( z i)) with 
a, a' = x, y. 



Losses 



The losses mainly due to the spontaneous emission 
from the atomic upper levels can be estimated by adding 
the imaginary parts as ~)l ss = 8HQ s /[n z (2\A ss \ + iT)] 
and ^ s - s = 4^[2 + cos(^-^)]n'/[n z (2A s5 + iT)], 
which leads to the total loss rate n = ^ s (k S s + k S s) 
with k ss = 8n ph , s n^r/ [n z (4A 2 ss +T 2 )] and k s - s = 

4 (2 + cos - <p a )) n ph , s H'r/ [n z (4A^_ + T 2 )} . With 
the smallest detunings A ss / ~ 4r used in the plotting 
of Fig. 2 for example, the shortest coherence time 1/k is 
400/is. For larger detunings, the coherence time will be 
longer. 
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